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soft supersymmetry-breaking couplings in a general softly broken N = 1 super- 
symmetric model. We also specialize these results to the Minimal Supersymmetric 
Standard Model. 



I. Introduction 



In the Standard Model, the mass of the Higgs scalar boson is subject to quadratically 
divergent radiative corrections which ought to be of order 5m? H ~ A 2 , where A is an ul- 
traviolet cutoff scale. This presents a naturalness problem if A 2 3> fnjj, since fine tuning 
is then necessary to explain why the Higgs mass is near the electroweak scale. Low-energy 
supersymmetry[l] evades this naturalness problem because the quadratic divergences can- 
cel. To accomplish this, the Minimal Supersymmetric Standard Model (MSSM) introduces 
a "sparticle" partner for each standard model particle. If supersymmetry were exact each 
sparticle would be degenerate in mass with its standard model partner, which is certainly 
not the case experimentally, so supersymmetry must be broken. Fortunately, the cancel- 
lation of quadratic divergences still works if supersymmetry is broken softly [2] by terms of 
dimension 2 and 3 in the Lagrangian. The masses of the sparticles are then determined 
by the soft super symmetry- breaking terms. 

Because low-energy supersymmetry is a perturbative solution to the naturalness prob- 
lem, one can attempt to relate observed phenomena at low energies to physics at very high 
energy scales. For example, it is remarkable that in the MSSM, the three gauge couplings 
appear to unify[3] at a scale 10 15 -10 16 GeV, hinting at a Grand Unified Theory (GUT) or 
some other organizing principle such as string theory. With the eventual discovery of the 
sparticles and determination of their masses, we should gain information about the soft 
supersymmetry-breaking parameters. Already we know that these parameters are not at 
all arbitrary, because otherwise large flavor-changing neutral currents would arise in the 
low-energy physics due to the effects of loops containing squarks, and the arbitrary com- 
plex phases of the soft supersymmetry-breaking terms would be expected to give rise to a 
CP-violating electric dipole moment for the neutron in violation of experimental bounds. 
Thus there is strong circumstantial evidence in favor of some organizing principle governing 
the soft supersymmetry-breaking terms. 

Models obtained from supergravity[4] can provide just such an organizing principle 
for the soft supersymmetry-breaking terms specified at some very high input scale. In 
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minimal supergravity, the spinless particles in the theory all obtain a common mass mo 
at this input scale. In this scenario, the absence of large flavor- changing neutral currents 
can be ascribed to the consequent near-degeneracy of the squarks. There are also scalar 
trilinear couplings among the squarks, sleptons and Higgs scalars as allowed by .R-parity 
and gauge invariance; these are each equal at the input scale to the corresponding Yukawa 
coupling times a universal mass parameter A. There is also a supersymmetry-breaking 
scalar (mass) 2 in the Higgs sector. Finally, there are three gaugino masses which also 
break supersymmetry. Large CP violation can be avoided with a common complex phase 
for all of the soft supersymmetry-breaking parameters. If one makes further assumptions 
about the high energy physics, additional constraints on the parameters of the theory are 
obtained. For example, if there is a GUT, then the gauge couplings and gaugino masses 
are unified, and there may be relations among the Yukawa couplings and among the scalar 
trilinear interactions at the unification scale. However, since the supersymmetry-breaking 
mechanism remains mysterious, a precise formulation of the organizing principle behind 
the soft terms in the MSSM remains unclear. 

In any case, with the parameters of the model specified at the input scale by some 
candidate organizing principle, one can run the couplings and masses down to low energies 
using the renormalization group equations and make predictions about the sparticle masses 
and other low-energy phenomena. Many authors (for example, [5-17]) have provided nu- 
merical and analytical results for the sparticle spectrum under a variety of assumptions 
and constraints on the input parameters. In studies of this type, greater precision as well 
as an estimate of errors incurred in the running can be obtained by employing the two- 
loop renormalization group equations. The two-loop /3-functions for the supersymmetric 
couplings (gauge couplings[18] and superpotential parameters [19]) have been known for 
some time. Recently, the two-loop /3-functions for gaugino mass parameters have also 
been found[20,21]. In this paper, we will complete the list of two-loop /3-functions for a 
general softly-broken supersymmetric model, by computing the results for scalar interac- 
tions which break supersymmetry softly. We will also give the results of specializing these 
calculations to the MSSM. 
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We consider a general N = 1 supersymmetric Yang-Mills model. The chiral superfields 
&i contain a complex scalar fa and a two-component fermion ^ which transform as a 
(possibly reducible) representation R of the gauge group G. The superpotential is 

W=±Y i i%Q j Q k + ±fJ%Q j + L% . (1.1) 

In addition, the Lagrangian contains soft supersymmetry-breaking terms of the form 

£SB = —h^fa^k - " l(m 2 )U% - \m\X + h.c. (1.2) 

where M is the mass of the gaugino A. The renormalization group equations for the gauge 
coupling and the superpotential parameters Y^ k , \s}i and U and the gaugino mass M are 
known. In this paper we give the corresponding results for the soft-breaking parameters 
W k , b l i, and [vr?)\. 

For simplicity we first give our results for the special case of a simple [or U(l)] gauge 
group. We then explain the modifications required if the gauge group is a direct product 
in Section III, and discuss the specialization of these results to the MSSM in Section IV. 
We let t A = (t)f^ denote the representation matrices for the gauge group G. Then 

(tV)j = C(R)5J (1.3) 
Tr^OA 5 ) = S(R)5 AB (1.4) 

define the quadratic Casimir invariant C(R) and the Dynkin index S(R) for the repre- 
sentation R. For the adjoint representation [of dimension denoted by d(G)], C(G)5 AB = 
f ACD f BCD with f ABC the structure constants of the group. 

In principle, the two-loop /9-functions for a general renormalizable theory have already 
been given in [22,23,24]. However, there are two practical issues which must be addressed 
in order to apply those results to the case at hand. 

The first issue is that only dimensionless couplings appeared in [22-24]. The two-loop 
/3-function for a gauge coupling was given in [22], for a Yukawa coupling in [23], and 
for a scalar quartic coupling in [24]. A general renormalizable theory may also contain 
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fermion mass terms, scalar (mass) 2 terms, and trilinear scalar couplings. Fortunately, the 
/3-functions for each of these dimensionful couplings can be inferred from the results given 
in [23] and [24] by taking some of the scalar fields to be non-propagating "dummy" fields 
with no gauge interactions. For example, a fermion mass term has the form 

C M = -M lj lpilpj = -^dummy^dummy^j • (1-5) 

Now if Vdummy i s taken to have no other interactions, then the /9-function for the fermion 
mass has the same form as that of the Yukawa coupling ^dummy' Thus the (3- function 
for any fermion mass can be inferred directly from the results of [23]. Similarly, scalar 
(mass) 2 terms can always be thought of as scalar quartic interactions involving two dummy 
scalars and two normal scalars, while scalar trilinear couplings can be thought of as scalar 
quartic interactions with one dummy scalar and three normal scalars, so the two-loop f3- 
functions for those cases can be inferred from the results of [24] using a small amount of 
careful combinatorics associated with the symmetry factors. 

The second issue to be addressed is that the results of [22-24] were obtained using 
dimensional regularization [25] (DREG). Now, DREG violates supersymmetry explicitly 
because it introduces a mismatch between the numbers of gauge boson and gaugino degrees 
of freedom. Therefore the use of DREG is inappropriate for (even softly broken) super- 
symmetric models. Instead, one should use the modified scheme known as dimensional 
reduction [26] (DRED) which does not violate supersymmetry. 

In DREG, supersymmetry is violated in the finite parts of one-loop graphs, and in the 
divergent parts of two-loop graphs. This means that for a given set of physical quantities 
(e.g. pole masses and S-matrix elements) the running couplings computed in DREG with 
modified minimal subtraction[27] (MS) will differ from those computed in DRED with 
modified minimal subtraction (DR) by finite one-loop corrections, and the /3-functions will 
be different for the two schemes starting at the two-loop level. We therefore present our 
/3-functions as they appear in the DR scheme. This means that the results of [22-24] must 
be translated from MS to DR. In [20] , we provided a "dictionary" for translating couplings 
between the two schemes including all finite one-loop radiative corrections, which is the 
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order necessary for the present application to two-loop /^-functions. (The relationship 
between DRED and DREG for non-supersymmetric theories has recently been illuminated 
in [28].) We must also take into account the fact that in DRED, the so-called e-scalars[26] 
obtain one-loop mass counterterms due to the supersymmetry-breaking scalar and gaugino 
masses. 

In summary, our method of computation is to use the results of [22-24] for a general 
renormalizable theory, using dummy scalars where necessary for dimensionful couplings, 
and translating the results from MS to DR. This process is straightforward but rather 
tedious, and we decline to exhibit the details here. It may also be possible in principle to 
extract these results from the dimensional reduction calculations of [29]. 



II. Two-Loop Running in a General Softly Broken Supersymmetric Model 

For completeness, we begin by reviewing the known results for gauge couplings, gaugino 
masses, and superpotential parameters (including Yukawa couplings). This will also serve 
as a useful point of reference for the notation established in the Introduction. We will 
then provide the two-loop beta functions for h^ k , , and (m 2 );-, which constitute the 
new results of this paper. We do not assume anything about the relative complex phases 
of any of the parameters. 

The gauge coupling at two loops is actually scheme-independent. Therefore, it may 
be obtained simply by specializing the results of [22] to a general supersymmetric model. 
Doing so, we obtain the known[18] result: 

-Q ~ 1 flW I 1 3 {2) (2 1) 

dt 9 ~ 16n^ g + (16tt 2 ) 2 ^ { } 

4 1} = 9 3 [S(R) - 3C(G)} (2.2) 
4 2) = g 5 {-6[C(G)] 2 + 2C(G)S(R)+4S(R)C(R)}-g 3 Y^% k C(k)/d(G) . (2.3) 

Here = (Y 1 -^)* and S(R) is the Dynkin index summed over all chiral multiplets and 
S(R)C(R) is the sum of the Dynkin indices weighted by the quadratic Casimir invariant. 
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The two- loop (3-f unction for the gaugino mass parameter may be computed (as we did 



in [20]) by first using the results of [23] for a general theory in MS and then translating 



the results to the DR scheme appropriate for supersymmetry. The result is: 



1m 

dt 



1 



1 



T0 



(2) 



16yr 2^M t- ( 167r 2)2^M 

0$ =g 2 [2S(R) - 6C(G)} M 

0$ = g 4 |-24[C(G)] 2 + 8C(G)S(R) + 16S(R)C(R)} M 



(2.4) 
(2.5) 



+ 2g< 



h ijk _ M yyAj Y ijk C(k)/d(G) . (2.6) 



This result was also obtained by Y. Yamada[21] using a different method. Note that we 
have adopted a slightly different notation here than we did in [20], since the heights of 
indices on the last term in (2.6) have been reversed. This is the same as exchanging the 
role of the gaugino mass M and its conjugate . We do this to ensure consistency with 
results below while maintaining the fewest possible number of M^s. 

The two-loop /3-functions for the superpotential parameters can be obtained either by 
superfield techniques[19], or by applying the general results of [23] and [24]. In applying 



the latter method, one must be careful to convert the MS results into DR, as we have 
already mentioned. We obtain the same results using both methods: 



dt 



d_ 

dt 
d 



ip 



2 IP 



1 



(16tt 2 ) 



2^2 > 



(2)k 



16tt 



2\2 7 P 



(2)j 



+ (k <-> i) + (k <-> j) 
+ U ^ i) 



—L l = LP 
dt 



2^ + 



167T 



(16tT 2 ) 



2\2 7 P 



(2)i 



where 



7. 



lf )] = -\Y imn Y nm Y pqr Y™3 + g 2 Yipq Y^[2C(p) - C(i)} 
+ 28{g 4 [C(i)S(R) + 2C(i) 2 - 3C(G)C(i)} 



(2.7) 
(2.8) 
(2.9) 

(2.10) 
(2.11) 



In these equations, C(r) always refers to the quadratic Casimir invariant of the representa- 
tion carried by the indicated chiral superfield, while S(R) refers to the total Dynkin index 
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(D? (2)?' 

summed over all of the chiral superfields. The objects 7' ; and 7. arise completely 
from wave-function renormalization in the superfield approach, in accordance with the 
"non-renormalization" theorems of supersymmetry [30] . 

Next we consider the two-loop /^-function for the soft supersymmetry-breaking scalar 
trilinear coupling h^ k . This is obtained by specializing eqs. (4.3) and (4.7) of [24] to the 
supersymmetric case, with one of the external scalar fields being a dummy field, and then 



translating the result from MS to DR using [20]. The result is: 



dt 



ijk 



16tt 2 



4 1 ' 



ijk 



+ 



ijk 



(16tT 2 ) 2 

Imn* f Y^ l Y lmn h mnk - 2 [h^ k - 2MY^ k ) g 2 C{k) 
+ (k <-> i) + (k <-> j) 



(2.12) 



(2.13) 



0(2) 



hijlv, V n P1V vmrk 
2 I lmn I tpqr* 

vijW, V n PQV h mr k vijlv, h n PQV vmrk 
1 ' I lmn I I pqr"' 1 1 lmn n 1 pqr 1 

+ (h^ l Y lpq YPl k + 2Y^ l Y lpq h pqk - 2MY i i l Y lpq Y pqk } g 2 [2C(p) - C(fc)] 
+ (2h ljk - 8MY 1 ^ g 4 [c(k)S(R) + 2C(k) 2 - 3C(G)C(k) 



+ (k <-> i) + (k <-> j) 



(2.14) 



Next we consider the two- loop /3-function for the scalar (mass) 2 This is again 
obtained by specializing eqs. (4.3) and (4.7) of [24] to the supersymmetric case, but with 
two of the external scalar fields as dummy fields, and then translating the result from MS 



to DR using the results of [20]. Doing so, we find: 



d 
dt 



b*3 



A 1} 



13 



+ 



(16tt 2 ) 



2^2 



4 2) 



(2.15) 



= \b ll Yi mn Y mn i + \Y^ l Y lmn b mn + » il Y lmn h mn j - 2 ( 6« - 2M^ ) g*C(i) 



+ j) 



(2.16) 



0, 



(2) 



- ^ l YimnY pqn Y pqr Y™l - -Y^ l Y lmn b mr Y pqr YP qn 



_ IaA^V. n mr V h pqn - n il Yi h npq Y V mr i 
2 1 lmnh b I pqr"' P 1 lmn IL 1 pqr 1 

- » il Y lmn Y npq Y pqr h mr j + 2Y^ l Y lpq (6M _ <7 2 C(p) 
+ (&* W^' + 2vPY lpq hPti - rt w W) # 2 [2C(p) - C(i)] 
+ (2&y - 8^m) / [c(i)£(i2) + 2C(z) 2 - 3C(G)C(z)" 
+ (i <-> j) . 



(2.17) 



Finally, we consider the two-loop /3-function for the scalar (mass)^ (m 2 ) J i . This is once 
again obtained by specializing eqs. (4.3) and (4.7) of [24] to the supersymmetric case, with 



two of the external scalar fields as dummy fields, and translating from MS to DR using 
[20]. One must also be careful here to take into account the mass counterterms for the 
e-scalars which arise in DRED. The result is: 



Jt [m 



111 



J 

J i 



1 



171 



+ 



1 



16tt 2 

\Y im Y^ n {m 2 ) 3 n + '-Y^nim 2 ) 7 ! + 2Y ipq Y^{m 2 ) q 



(167T 2 ) 2 
1, 

2 

+ K pq h m - 8S^MM^g 2 C(i) + 2g 2 tf J Tr[t A m 2 ] 



(2.18) 



(2.19) 



P {2 1 

771 



' = -^(m 2 )^ ?mn y^y rer y^-i(m 2 )/y^y mrj y^y p(zn 

l z z 

- Y llm Y^(m 2 ) l r Y npq Y^ - Y Um Y^(m 2 ) r n Y rpq Y^ 

_ Yii m Y^ nr (m 2 ) n Y pqr Y pqm — 2Yii m Y^ n Y npq Y mpr (m 2 ) r 



(m 2 )k Pq y JPq + Y ipq Y lpq (m 2 )j + AY ipq Y^ l (m 2 ) q l + 2h ipq W pq 



g 2 [C(p) + C(q) - C(i)} 



~~ Yn m Y^ n h npq h mpq — hii m h^ n Y npq Y mpq 

~ hu m Y^ n Y npq h mpq — Yn m h^ n h npq Y mpq 

r. o.z__ 

+ 

- 2h ipq Yi pq M - 2Y ipq h? pq M^ + AY ipq Yi pq MM^ 
- 2g 2 tf\t A m 2 ) l r Y lpq Y rpq + 8gh? j Tr[t A C(r)m 2 } 
+ 5{g A MM^ \2AC{i)S{R) + 48C(z) 2 - 72C(G)C{i) 
+ 8Sig 4 C(i)(Tr[S(r)m 2 } - C(G)MM^) . (2.20) 

Here h^y. = (h l]k )* . The traces are over all of the chiral superfields, and the C(r) are the 
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quadratic Casimir invariants for the irreducible representations of chiral superfields in the 



This completes the list of two-loop /3-functions for a general softly-broken supersym- 
metric theory. We conclude this section by noting a non-trivial consistency check on these 
results. It has been shown[31] that an N = 2 supersymmetric Yang-Mills theory is finite 
to all orders in perturbation theory provided that certain constraints are imposed on the 
representations of superfields. This finiteness continues to hold even with soft breaking, 
provided that the soft terms obey certain additional constraints [32]. It is an amusing ex- 
ercise to check that when the above formulas for a general N = 1 supersymmetric model 
are specialized to the N = 2 case with the appropriate contraints on the representations 
and soft couplings, the two-loop /9-functions do indeed vanish. 



As promised, we now point out the modifications which must be made to the preceding 
formulas if the gauge group is a product of simple [or U(l)] subgroups G a - 



One obtains the /3-function for each gauge coupling g a by applying the following rules 
to (2.2)-(2.3): 



traces. Note that the terms which explicitly involve t i 



are zero for non-abelian groups. 



III. Two-Loop Beta Functions for Direct Product Groups 



g 3 C(G) - glC{Ga) 
g' 3 S(R) -> g 3 a S a (R) 
g 5 C(G) 2 - glC(G a f 
g 5 C(G)S(R) - glC{G a )S a {R) 



(3.1) 



(3.2) 



(3.3) 



(3.4) 




(3.5) 



b 



g 3 C(k)/d(G)^glC a (k)/d(G a ) . 



(3.6) 



The J2b i n (3-5) is a sum over subgroups. Similarly, one obtains the /3-function for each 
gaugino mass parameter M a by applying to (2.5)-(2.6) the same rules as above, but with 
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one less power of g — > <7 a , and with M — > M a wherever it appears. The exception is that 

16g 4 S(R)C(R)M ^8j2d 2 a9b S a( R ) C b(R) (M a + M b ) (3.7) 

b 

in (2.6). Note that in these cases, all terms which do not contain a quadratic Casimir 
invariant of a non-adjoint representation are simply diagonal in each subgroup. 

For the (3- functions of the superpotential parameters and soft supersymmetry-breaking 
scalar interactions (2.7)-(2.20), one always obtains a sum over subgroups. In each of these 
cases, the following set of rules apply: 

9 2 C{r) ^Y,9lCa{r) (3.8) 

a 

g 4 C(r)S(R) -> £ giC a (r)S a (R) (3.9) 

a 

g A C(r)C(G) -> £ ^C a (r)C(G a ) (3.10) 

a 
a 

Terms which also involve gaugino masses are modified exactly as above, with 

M, Aft — > M , Af] (3.12) 

except for the term 

48/MMtc(z) 2 ^ X! E [32M a M a f + 8M a M 6 f + 8M 6 mJ] (3.13) 

a 

in eq. (2.20). We also have in (2.19)-(2.20): 

gh*Tr[t A m 2 ] - £ <7 2 (tft'Tr[t fm 2 ] (3.14) 
a 

a 

/t^Tr[t^C(r)m 2 ] - £ £ A 6 W)^[t£ Q,(r)m 2 ] (3.16) 

a 

/C(z)Tr[S(r)m 2 ] - £ ^C a (z)Tr[S a (r)m 2 ] . (3.17) 

a 
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Finally, we must note that there is an exceptional case when the gauge group contains 
(unlike the MSSM) a direct product of more than one U(l). One should then choose the 
basis for the U(l) subgroups so that the matrix Tt[q a qfj\ is diagonal, where the trace is 
over all chiral superfields and q a denotes the U(l) a charge. This is always possible, since 
Tr[(/ a (/5] is always a real symmetric matrix. Then the only non-trivial rule is that the 
term g 4 C(i)Tr[S(r)m 2 } in (2.20) becomes a sum over non-U (1) subgroups as before, plus 
a contribution 

/C(z)Tr[S(r)m 2 ] -> £ £ g&g(q a ) M) ;Tr[ M m 2 ] (3.18) 

a b 

where J2 a and are sums over U(l) subgroups, and (q a )i denotes the U(l) a charge of 
the chiral superfield carrying the index i. 

IV. Two-Loop Running in the Minimal Supersymmetric Standard Model 

In the MSSM, the gauge group is SU(S) C x SU(2) L x U(1) Y , with chiral superfields Q 
and L for the 5'L r (2) L -doublet quarks and leptons, and u, d, e for the 5L r (2) i -singlet quarks 
and leptons, and two Higgs doublet chiral superfields H u and H d . The superpotential is 

W = uY u QH u + dY d QH d + eY e LH d + fiH u H d (4.1) 

where Y u , Y d , Y e are each 3x3 Yukawa matrices. The soft supersymmetry-breaking 
Lagrangian contains scalar couplings 

-C = iih u QH u + dh d QH d + eh e LH d + BH u H d + h.c. (4.2) 

where h u , h^, h e are again 3x3 matrices in family space, and a hat is used to denote the 
scalar component of each chiral superfield. There are also scalar masses of the (m 2 )^ type: 

-C = m 2 Hu HlH u + m 2 Hd H d H d + Q^ytiqQ + L^m 2 L L + um^v) + dm 2 d S + em^ . (4.3) 

Here again mg, m 2 , m 2 , m 2 , and m 2 are 3x3 matrices in family space. Finally, the 
gauginos for the subgroups SU(3) C , SU(2) L and U(1) Y have masses M3, M2, and Mi, 
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respectively. Again, we do not assume anything about the complex phases of any of the 
parameters. For most applications, it will be sufficient to retain only the Yukawa couplings 
of the heaviest family, but we prefer to retain complete generality. 

For the sake of completeness and to provide a useful point of reference, we begin by 
reviewing the known results for two-loop /^-functions within the MSSM. For the three 
gauge couplings, we have from (2.1)-(2.3): 



o 3 



dt 9a ~lto? a + (16tt2)2 



6=1 



x~u,d,e 



(4.4) 



,(1) 



Here B d ' = (33/5,1,-3) for U(1) Y (in a GUT normalization), SU(2) L , and SU(3) C re- 



spectively, and 



^ = 1 »/? 



/ 199/25 27/5 88/5 \ 



25 
9 



24 

14 J 



and 



a 



u,d,e 



'26/5 14/5 18/5' 
6 6 2 
4 4 



(4.5) 



V H/5 

The two-loop renormalization group equations for the three gaugino mass parameters 
[20,21] can then be written easily in terms of the same coefficients: 

V ' 6=1 ^ (4.6) 

+ C^(Tt[Yhx]-M a Tt[Y^Y x ^ 

x=u 1 d 1 e 



From (2.7)-(2.11), the two-loop beta functions for the superpotential parameters are: 



d 
dt 



d 



167T 
1 



— + 



(IGtt 2 ) 2 ^ 



(2) 



+ 



r/5- 



(2) 



with: 



dt u ^ e IQ^u,^ ^ ( 167r 2 )2 ^Y„,^ e 
/#> = ^{Tr(3Y u Yl + 3Y d Y+ + Y e Yj) - 3g 2 2 - 



(4.7) 
(4.8) 

(4.9) 



fi{ - 3Tr(3Y u YlY u Yl + 3Y d Y+Y d Y+ + 2Y U Y^Y^Y^ + Y e YjY e Yj) 
+ [l6^ 2 + ^]lt(Y tt Yt) + [l6^ 3 2 - ^ 2 ]Tr(Y d Yj) + ^ 2 Tr(Y e Y|) 



15 4 9 2 2 207 4 i 
+ Y 92 + 19192 + -^91 1 



(4.10) 
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$1 = Y u {3Tr(Y u Yl) + 3YiY u + Y\Y d - ^g 2 - Sg 2 - ^g 2 } 
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(4.11) 



A it 



Y u { - 3Tr(3Y u YlY u Yi + Y u Y\Y d Yl) ~ YjY d Tr(3Y d Y+ + Y 6 y|) 

- 9YiY u Tv(Y u Yi) - 4YiY u YlY u - 2Y^Y^Y^Y^ - 2Y^Y^Y^Y t 



+ 



i%3 + pi 



Tr(Y M Yl) + 



6^1 + \q\ 



16 4 o 2 136 2 2 15 4 2 2 
y#3 + 8 #3#2 + ~^9z9\ + y# 2 + 9 2 9i + ^ 



2743 4 



YlY M + Y^Y^ 



Si 



(4.12) 



(3^ d = Y d {Tr(3Y rf Yj + Y e Yl) + 3YjY (i + YlY 1 



fal-^-^al} (4-13) 



(2) _ 



d{ - 3TY(3Y d YjY d Yj + Y^Y^Y^Y^ + Y e YjY e Y|) 

- 3YlY u Tr(Y u Yl) - 3Y^Y^Tr(3Y^Y^ + Y 6 yJ) - 4Y^Y^Y^Y^ 

- 2YlY u YiY u - 2YiY u YtY d + \l6g 2 - ^g 2 } TV(Y d Y+) 



a 2 , ^ 2 
%2 + 7f 1 



+ ^ 2 Tr(Y e Yl) + ^ 2 YiY u + 

16 4 n22 §22 15 4 92 287 4I 
- + 8 ^2 + g^3^1 + y#2 + ^2^1 + -gQ-01 ) 



(4.14) 



= Y e {Tr(3Y d Yj + Y c yJ) + 3yJy 6 - Zg\ - ®-g 2 } (4.15) 

= Y e { - 3TY(3Y d YjY d Yj + Y u Y+Y d Yl + Y e Y|Y e Y|) 

- 3Y!Y e TY(3Y d Yj + Y 6 yJ) - 4YlY e YjY e + [l6</§ - Tr(Y^Y^) 

+ ^ 2 Tr(Y e Yt) + 6, 2 YjY e + ^g\ + \g\g\ + ^gf} (4.16) 



Of course, the /3-functions for Y u d e are 3x3 matrices in family space. 
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The above results for the MSSM have all appeared before. Now we apply our results 
of Sections II and III to arrive at the two-loop beta functions for the soft-breaking trilinear 
scalar couplings: 

rl 1 n\ 1 

(4.17) 



d 1 (i) 



1 -3 {2) 



(3^ = h u {3Tr(Y u Yj i ) + 5YiY u + YjY rf - ^-g 2 - 3g 2 - ^g 2 ) 

+Y M {6Tr(h u Yl) + AYih u + 2Y\h d + y^M 3 + Qg 2 2 M 2 + ^<7 2 Ml} (4.18) 



13 



(3 { H = h u { - 3Tr(3Y u Y ? 1 i Y u Yl + Y.yJy^yI) - YjY d Tr(3Y rf Yj + Y e Y|) 

- 15YlY u Tr(Y u Yi) - 6y1y u yIy u - 2Y^Y^Y^Y^ - 4Y+ Y^Y, 



+ 



16^ + ^g\] Tr(Y u Yl) + 12^ 2 y1y u + |^Y+Y d 
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Lt) 4 . o 2 2 . 136 2 2 , ib 4 . 2 2, ZY4a 4\ 
+ 8 #3#2 + ~^-9s9l + y 92 + 929l + ~^q9i j 

+Y U [ - 6Tr(6h u YlY u Yl + h^Y^Y* + h^Y^Y*) 

- 18YlY u Tr(h M Yl) - Y+ Y d Tr(6h d Yj + 2h e Yj) - 12Ylh u Tr(Y u Yl) 

- Yjh d Tr(6Y d Yj + 2Y e Y|) - GY^Y^ - 8Yyh u YyY M 

- 4YjY d Yjh d - 4Yjh d Yt Y d - 2Y^ Y^yJ h u - 4Y> d Yt Y u 
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2743 



+ 



32s| + Tr(h u Yy) + \&g 2 2 + \g 2 ] Y+h u + ^?Yjh d 



32^M 3 + -^ 2 Mi]Tr(Y u Yb 
64 



12g%M 2 + \g\M x 
272 



Y 1 Y 

1 11. 1 



u 1 u 



-g^M^Ya + -glM 3 - 16<7 3 2 <7 2 2 (M 3 + M 2 ) - — <7 3 V(M 3 + M x ) 



- 30glM 2 - 2g 2 g 2 (M 2 + M x ) - ^9^] 



(4.19) 



$ = h rf {Tr(3Y rf Yj + Y e Yj) + 5YjY d + yIy u - -g 2 - 3g 2 - -g 2 ] (4.20) 
+Y rf {Tr(6h d Yj + 2h e Yl) + 4Yjh d + 2Yjh u + fg$M 3 + 6g 2 M 2 + ^<7 2 Ml} 

tf® = h d { - 3Tr(3Y rf YjY rf Yj + Y u Y rf Y d YU Y e Y|Y e Y|) 

15 



3YiY u TT(Y u Yl) - 5Yj Y d Tr(3Y d Yj + Y e Yj) - GY^Y^ 



2Yy Y u y}jY u 



4YlY u YtY d + \l6gl - -^2]iv(Y rf Yj) 



+ §<7 2 Tr(Y e Y|) + ^?YtY u + [l2^ 2 2 + ^g 2 ] Y t Yrf 



5 l 

16 4 
.- ■ ' 



2 2 8 o o 15 4 2 2 287 4~| 
- 1T^3 + %3#2 + g939l + y #2 + 9 2 9l + y^l / 



90 



+Y d { - 6Tr(6h d YjY rf Yj + h^Y^Yl + W^Y^ + 2h e Y|Y e Y|) 

- 6Y] 1 Y u TV(h u Yj i ) - 6YjY rf Tr(3h d Yj + h e Y+) 

- 6Ylh u IY(Y u Yi) - 4Y> d Tr(3Y d Yj + Y e Y|) - 6YjY rf Y> d 

- 8Ylh,Y),Y d - 4Ylh w YlY u - 4Yj i Y u Yj i h u - 4Yjh u YlY 



2YlY M Yjh rf + [32g 2 3 - Tr(h d Y+) + ^^ 2 Tr(h e Yl) + j^Y^ 



+ 



%2 + 



32^M 3 - -^Mill^Y^) - — ^MiTr(Y e Yt) 



12^M 2 + -g\M x 



Y d Y d 



IqImMYu + y^M 3 - 16^ 2 2 (M 3 + M 2 ) 



y <7§<7?(M 3 + Mi) - 30gjM 2 - 2g 2 2 g 2 (M 2 + M x ) - ^ 



01 Mi} 



(4.21) 



^ = h e {lY(3Y d Y+ + Y e Ye) + 5yJy 6 - 3g 2 2 - jj<? 2 } 

+Y e {Tr(6h d Y^ + 2h e Yj) + 4Yjh e + Qg 2 2 M 2 + y £ 2 Mi} 



(4.22) 



= h e { - 3Tr(3Y d YtY d Yj + Y^Y^Y^Y^ + Y e YjY e Y|) 

- 5YjY e TY(3Y d Yt + Y e Yj) - 6Y|Y e YjjY e + [l6^ 3 2 - 2 -g 2 } Tr(Y^Yt) 



+ ^ 2 IV(Y e Yt) + [l2, 2 - \g 2 ] YtY e + ^g\ + \g\g\ + ^gf} 
+Y e { - 6Tr(6h d YjY rf Yj + KY^Y^l + h^Y^Y* + 2h e YjY e Yj) 
- 4Yjh e IV(3Y d Yj + Y e Y|) - 6YjY e Tr(3h d Yj + h e Y|) 



- 6YjY e Yih e - 8Ylh e YjY e 
4 o1„ ,. 12 



+ 



32<7 3 2 - -g\ TV(h d Yj) + -^ 2 IY(h e Y|) + 6^ 2 + -a 2 



Yjh, 



32g 2 M 3 - \glMM^{Y d Y\) - ^<7 2 MiTr(Y e Y|) - 12^ 2 2 M 2 YjY e 
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30^M 2 - y(&J(Mi + M 2 ) - 54^} 



(4.23) 



These are again 3x3 matrices in family space. One should note that there is no universally 
agreed-upon convention for the sign of the gaugino masses in these equations, because of 
the freedom to rotate the phases of the gaugino mass terms in the Lagrangian. 

The MSSM also contains one scalar (mass) 2 of the type which from (2.15)-(2.17) 
satisfies the two-loop renormalization group equation: 



d 
~dt 



B 



+ 1 *(2) 



16tt 2 ^ t (16tt 2 ) 2 ^ 



(4.24) 



0$ = s{Tr(3Y u YU3Y rf Yj + Y e Y|) - 3<^ - |</f} 

+^{Tr(6h u Yi + 6h d Y+ + 2h e Y+) + 6^M 2 + ^ 2 Mi} 



(4.25) 



/4 2) = S{ - 3Tr(3Y u Yj i Y u Yl + 3Y^Y^Y^Y^ + 2Y U Y^Y^Y^ + Y e Y|Y e Y|) 
+ [l6^ 3 2 + ^ 2 ]Tr(Y u Yl) + [l6^ 3 2 - ^Tr^Y*) + ^ 2 Tr(Y e Yj) 



15 4 9 2 2 207 4I 
+ Y 92 + 5^2 + ^1 1 



+^{-12Tr(3h u YlY u Y] i + 3h rf YjY rf Yj + h u Y+Y d Yl + h^Y^Y* + h e YjY e Yl) 



+ 



32^ + -g\] Tr(h u Yl ) + [32«/§ - -g\] Tr(h d Y+) + -^ 2 Tr(h e Yj) 
32^ 2 M 3 + ^M^TviYuYt) - \32g 2 M 3 - I^mJ Tr(Y d Yj) 



- y^MiTr(Y e Y+) - 30^M 2 - |4 2 (Mi + M 2 ) - ^<7l^l} • (4-26) 
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Finally, we turn to the /^-functions for the scalar (mass) 2 terms of the (m 2 )j type in 
the MSSM. It is convenient to define the quantities 



S = m H u ~ m H d + Tr t m Q - m i - 2m l + m d + m e 



(4.27) 



and 



S' =Tr 



iSm 2 Hu + m 2 Q )YlY u + 4Ylm 2 Y u + (Sm^ - m^Y+Yd - 2Y rf m 2 Y d 
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and 



+ 



+ (m 2 Hd + m 2 L ) YlY e - 2Yjm 2 Y e ] 
{™ 2 H u - m H d - Tr ( m l)} + 



3 o 3 o 
^2 + ^1 

16 2 16 2 
y^3 + 15^1 



Tr(m 2 ) + 



5 2 , z 2 
3^3 + 15^1 



3^3 + ^2 + 3^1 



Tr(m 2 ) + ^ 2 Tr(m 2 ) 



lV(m|) 



(4.28) 



0"2 
^3 



^?{3(m| u + + Tr[mg + 3m| + 8m 2 + 2m 2 + 6m 2 ]} 

92 {™ 2 H U + ™H d + TV[3m| + m|] } 
^Tr[2m| + m2 + m 2] . 



Then from the results of sections II and III, we obtain for each of the /^-functions in 
the standard form 

^m 2 = * ^\ + -^-J 2 ] (4.29) 

dt 167T 2 m (167T 2 ) 2 m V ' 

the following results: 



rnir 

tin 



= 6Tv[(m 2 Hu + m 2 Q )YlY u + Y^m 2 Y u + h+h, 

_ 6 ^2| M2 |2_6 2| Mij 2 3 2 S 



(4.30) 



' mir 

tt>U 



- 6Tr [6(m^ + m^) y!y u yIy u + GY^m, 2 Y u yIy u 

+ (™>H U + ™ 2 H d + m^)YiY tt YjY d + Ylm 2 Y ?i YjY d 
+ Y^m^Y* Y d + Y^Y^m 2 Y d + Gh^Y^ + G^Y^Y^ 
+ hJh rf YlY u + Y\ Y d hUi + hjY d Yth u + YjhrfhiY* 
32^ + Tr[(m^ u + m 2 Q ) Y^Y^ + Y^m 2 Y u + hlh u ] 
32^{2|M 3 | 2 Tr[YlY M ] - M 3 *Tr[Yih M ] - M 3 Tr[h|X]} 
^ 2 {2|M 1 | 2 Tr[YlY u ] - MfTr[Ylh M ] - MiTr^Y^} + ^g 2 S' 



+ 
+ 
+ 



+ 33^ 2 4 |M 2 | 2 + y^ 2 V(|M 2 | 2 + |M X | 2 + Re[MiM 2 *]) + \M X \ 2 



(4.31) 
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Tr 



6(m 2 Hd + m|)YjY d + 6Yjm 2 Y d + 2(m^ + m 2 )y|y 6 + 2Yjm 2 Y e 
+ 6h+h d + 2hjh e ] - 6s 2 2 |M 2 | 2 - ^i 2 |M!| 2 - (4.32) 



PI 



" 6Tr [6(r4 d + m|) yJy.yJy, + GY^m 2 Y rf Y, Y rf 

+ (™k + ™ff d + mg) YiY u YjY d + Y^m 2 Y u yJy^ 

+ YlY.m^YjYrf + YiY u Yjm 2 Y d + 2(m^ + m 2 )YjY e YjY e 

+ 2Yjm 2 Y e YjY e + Gh^Y* Y d + 6h+ Y rf Y*h d + h^Y* Y d 

+ YiY u hJh d + hlY u Yjh d + Y* h u h+ Y d + 2hJh e YjY e + 2hjY e Y|h e ] 



+ 
+ 

+ 



32^1 



#1 



n(m 2 Hd + m^YjY, + Ytm 2 Y rf + hjhrf] 



32^ 2 {2|M 3 | 2 IY[YtY d ] - M 3 *Tr[Y> d ] - M 3 Tr[h+Y rf ]} 
^ 1 2 {2|M 1 | 2 Tr[YjY d ] - MfTrYjhrf] - MiTr[h rf Y d ] | 
^ 2 {Tr[(m^ + m 2 )YjY e + Yjm 2 Y e + hjh e ] + 2|M!| 2 Tr[YjY e ] 
- M!Tr[hjY e ] - MfTYYjhe]} - ^ 2 <S' + 33^|M 2 | 2 



+ ^r9bl(\M2\ 2 + m\ 2 +Re[M 1 M 2 *]) + ^^|M!| 2 



+ 3# 2 a 2 + ^ 2 ai 



(4.33) 



The ,3-functions for nig, m 2 , m 2 , m^, m 2 are of course 3x3 matrices: 



P. 



(1) 

9 

ml 



(m 2 , + 2m% u )YiY u + (m 2 Q + 2m^)Yj Y d + [Y+Y, + Y^m 2 , + 2Ylm 2 Y t 
+ 2Yjm 2 Y d + 2hih u + 2hjh d - f <7 3 2 |M 3 | 2 - 6^ 2 |M 2 | 2 - ^-g\\M^ + \g 2 S 



(4.34) 



(2m 2 , + Sm 2 Hu )YlY u YlY u - 4Ylm 2 Y u YlY u - 4Yj f Y M mgY] i Y u 
4YlY u Ylm5Y u - 2Y] i Y u YlY u m 2 ) - (2m 2 , + 8m|j yJy^yJy^ 
4Ytm 2 Y dY t Y d - 4Yj Y d m£ Y* Y d - 4Y rf Y d Y^m 2 Y d - 2YjY d YjY d m 2 ? 
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"(m^ + 4m 2 H jYlY u + 2Ylui 2 u Y u + Y^m^ Tt(3yIy u ) 
~(m 2 Q + Am 2 Hd )Y\Y d + 2Yt m 2 Y d + y\ Y d m 2 Q ] Tr(3Yt Y d + yJy 6 ) 

- 6YlY u Tr(m^YlY u + Y^m 2 Y u ) 

- Y\ Y d TY(6m|Yt Y d + 6Y]p$Y d + 2m 2 yIy £ + 2Ytm 2 e Y e ) 
~~ 4|YyY u hJ[h u + h^h M YyY M + Y^huh^Y^ + h^YuY^hu | 
" 4 { Y d Y d^d + h> d Yt Yd + Yjh d ht Y d + htY.Yjh, 

- hih u Tr[6YiY u ] - Y* Y u TV[6hih u ] - h^Ti^Y^] - Ylh u Tr[6h+Y u ] 

- hJh d Tr[6Yt Y d + 2YjY e ] - Y^^h^ + 2hjh e ] 

- h d Y d Tr[6Yjh d + 2Yjh e ] - Yjh/fr[6ht Y d + 2hjY e ] 

+ ^{(2m2 + 4m^ u )YiY u + 4Ylm 2 Y u + 2YyY M mg + 4hlh M - 4MihiY u 
- 4MfYih u + 8|M!| 2 y1y u + (mj + 2m^)YjY d + 2Yjm 2 Y d 
+ Y d Y d m 2 Q + 2h d h d - 2M 1 h d Y d - 2AffY+li rf + 4|M!| 2 Y+Y d } 

+ ^ _ ^|| M3 |2 + 32 ^2 ( | M3 |2 + | Ma |2 + Re[ M 2 M 3 *]) 
32 

+ -^2 ()M3 |2 + | Mi) 2 + Re[MiM 3 *]) + 33^|M 2 | 2 
+ ^ 2 (|M 2 | 2 + |M X | 2 + He[MiM 2 *]) + ^7i|Mi| 2 

+ y^3°"3 + 3^2°"2 + hi (4-35) 



= (mf + 2mi )YjY e + 2Yjm 2 Y e + YjY e m£ + 2hlh 



6«i|M 2 | 2 -|</?|Mi| 2 -?</j5 



(4.36) 



/3 (2 1 = 



(2m 2 + 8m| d )Y|Y e Y|Y e - 4Y^m 2 Y e Y|Y e - 4YjY e m 2 yJy 6 



- 4y1y p y1iii 2 Y, - 2YjY P YjY P m 2 r 



(m 2 + Am 2 Hd )YlY e + 2Ylm 2 Y e + Y^m 2 ] Tr(3YjY d + y|y 6 ) 
YlYeTrtem^YjYrf + 6Yjm 2 Y d + 2m 2 Y+Y e + 2Y|m 2 Y e ] 

20 



- 4{YlY e h!;h e + hJh e YjY e + Y+hehjYe + hjY e Ylh e } 

- hJh e Tr[6YjY rf + 2YjY e ] - YjY e TV[6h^h d + 2h|h e ] 

- hjY e Tr[6Yjh d + 2Yjh e ] - Yjh e Tr[6h^Y d + 2hjY e ] 

+ ^! 2 {(m| + 2m| fd )Y e t Y e + 2Yjm 2 Y e + Y^m 2 + 2hjh e 

- 2M 1 h f e Y e - 2M 1 *Yjh c + 4|Mi| 2 Y|Y e } - ^-g\S' 
+ 33gl\M 2 \ 2 + ™g%gl{\M 2 \ 2 + \M X \ 2 + Re[MiM|]) + ^giW 
+ Sgfr2 + 501*1 (4.37) 



(2m 2 +4m|jY + 4Y u m 2 ? Yi + 2Y ?i Yim 2 i + 4h u hi 

-ffiSiMsi 2 -^^-^ 



(4.38) 



- (2m 2 + 8m|jY u YlY u Yi - 4Y u m 2 ? YlY u Y] 1 - 4Y u Ylm 2 Y u y1 

- 4Y u YlY u m 2 ? Y] i - 2Y u YlY u Ylm 2 - (2m 2 + 4m 2 Hu + 4m 2 Hd ) Y^yJy^yI 

- 4Y u m 2 ? YjY d Y] i - 4Y u Yjm 2 Y d Yl - 4Y„yJ Y d m£ Y* - 2Y u yJ Y^Y* m 2 

(m 2 + 4m^)Y u Yi + 2Y u m 2 ? Y] i + Y^m 2 ] Tt[6yIy u ] 

- 1 2 Y u Y^Tr [mg Y^ Y u + Y^m 2 Y u ] 

4 ^h u hJ[Y u Yjj + YyY^hyhy + h. u YyY u h]i + Y^hyhyY^ ^ 

- 4 |h u h)Y d Yi + Y u Yjh d hl + h u YjY d hl + Y^h^Y* } 

- 12 |h u hilV[YlY u ] + Y u YjlV[hih u ] + h u Yll¥[hiY u ] + Y u h£Tr [yJx] } 

+ [6^1 - 2 -g\] {(m 2 u + 2m 2 Hu )Y u Yl + 2Y u m 2 ? Yl + Y^m* + 2h u hl} 
12^{2|M 2 | 2 Y u Yl - M 2 *h M Yi - M 2 Y u hi} 

^?{2|M 1 | 2 Y ?i Yl - Mfh u Yi - MiY u hI} - ^g\S' 

(4.39) 



+ 



- ™ 9 lm\ 2 + 5 -§9hHm\ 2 + m\ 2 + Re[ Ml Mi}) + 3424 ^ »*- ' 2 



3 

16 2 16 2 



75 
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= (Ml + 4m^)Y d Yt + 4Y«,m^Yt + 2Y d Y*m5 + 4h rf h* 



/? (2 > = - 



f^l^3| 2 -^|M 1 | 2 + ^ 1 2 5 (4.40) 

(2m^ + 8m^)Y d YtY d Yt - 4Y d m 2 ? YtY (i Yi - 4Y^Y^m 2 Y^Y^ 
4Y,Yj Yrf m|Yt - 2Y d Yt Y.Ytm 2 - (2m 2 , + 4m% u + 4m 2 Hd )Y d Yl Y u y\ 

4Y rf m|YlY w Yj - 4Y d Yim 2 Y u y\ - 4Y d Yl Y u ui 2 q y\ - 2Y d YlY u Y d m 2 d 

i / i m ,2 \v .vt i ov ™2 vt i \/ _vt™2 rn,/ C vtv . i ivfv ~\ 



- - u — d u u—y, u d u u W u u " 

- ~(m 2 + 4ro^)Y d Yj + 2Y d m 2 ? Yt + Y^m 2 ] TY(6Y d Y d + 2Y|y, 

- 4Y d YjTr(3m|Yt Y d + SY^m 2 Y d + m 2 YjjY e + Ylm 2 Y e ) 

- 4{h,ht Y dY t + Y,Yjh d ht + h d Yt Y(i ht + Y^Y* } 

- 4{h d htY M Yj + YrfYlh^hJ + Y u hJ + Y^h^U 



4{h d htY u Yj + Y d Ylh u hJ + h d YiY u hJ + Y d hth u Yj} 
4h d htlY(3YjY d + YjY e ) - 4Y d YtTY(3hth d + hjh e ) 
4h d YtlY(3hjY d + h* Y e ) - 4Y d hJlY(3Yjh d + Yjh e ) 
%2 2 + §</?] { (™d + 2 m 2 Hd )Y d Yl + 2Y d m 2 ? Y* + Y^m 2 + 2h d hJ } 

• 12s 2 {2|M 2 | 2 Y d Yj - M 2 *h d Yj - M 2 Y d h^} 

- ^ 2 {2|M 1 | 2 Y d Yj - M*h d Yj - MxY d ht } + tfitf 

• ^ll^sl 2 + ^i 2 (|M 3 | 2 + \M X \ 2 + RefM^H) + 8 -^gj \M X \ 2 
16 o 4 9 

■ f7o (TQ -I f7i (Tl 



= (2m 2 + 4m 2 ) Y e Yj + 4Y e m 2 y| + 2Y e Yjm 2 + 4h e hJ 



-f^|M 1 | 2 + ^ 1 2 5 (4.42) 



^ = - (2m 2 + 8m 2 Hd )Y e YtY e Yt - 4Y e m 2 y!y 6 yJ - 4Y e Ylm 2 Y e Yl 
- 4Y e YjY e m 2 Y| - 2Y e YjY e Yjm 2 
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~(m 2 + 4m^)Y e Yj + 2Y e m 2 yJ + Y e Yjm 2 ] Tr[6Y+Y d + 2Y|Y e ] 

- 4Y e YjjTr [3mg Y^Y^ + SYjm^Y^ + m| yJy 6 + Y^m^] 

- 4{h e hjY e Yj + Y e Yjh e hJ + h e Y|Y e hJ + Y e hJh e Yj} 

- 4h e hjTr[3YjY d + Y+Y e ] - 4Y e Y|lY[3hJh d + hjh e ] 

- 4h e YlTr[3h^Y d + hjY e ] - 4Y e hjTr[3Yjh d + Yjh e ] 

+ [6^ 2 2 - ^i 2 ] {(m 2 + 2m 2 Hd )Y e Yt + 2Y e m 2 y| + Y e Yjm 2 + 2h e hJ} 
+ 12<7 2 2 {2|M 2 | 2 Y e Y3; - M 2 *h e Yl - M 2 Y e hJ} 

- ^ 2 {2|M!| 2 Y e Y| - Mfh e Yl - MiYehJ} 



o _/ 2808 4, , , ,9 12 9 



(4.43) 



V. Conclusion 

In this paper, we have presented the two-loop renormalization group equations for 
all couplings in a general softly-broken supersymmetric model, and in particular for the 
MSSM. If the sparticles predicted by the MSSM are found and their spectrum is determined 
with some accuracy (for example at an e + e~ collider [33]), these results may be useful in 
discriminating between various candidate organizing principles for the soft supersymmetry- 
breaking terms at some very high input scale. One can run the parameters down to 
low energies to predict the masses of the sparticle spectrum and other features of low 
energy phenomenology in terms of what may turn out to be only a few input parameters. 
The masses of the sparticles depend primarily on just those soft supersymmetry-breaking 
couplings whose /3-functions have been given to two loops here. We find that the two- 
loop /3-functions generally make a difference of several percent (compared to the one-loop 
predictions) for the sparticle masses, although it is quite difficult to make meaningful 
estimates of the size of the two-loop corrections without committing to a specific model. 
At the same level of accuracy, one must also be careful to take into account threshold 
effects as well as the distinction between running masses and pole masses. In extensions of 
the MSSM which have a large non-minimal particle content above the electroweak scale, 
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the two-loop corrections and threshold effects are potentially much larger. 

Note added. Since the original preprint version of this paper appeared, the (3 functions 
of Sec. II have been calculated by Yamada [34] using superfield techniques, and by Jack 
and Jones [35] working directly in DRED. Our results agreed with theirs except for dis- 
crepancies in the two-loop (5 function (2.20) for scalar masses of the (m 2 )] type, namely 
the coefficient of the last term proportional to 5lg A C{i){Tx[S{r)m 2 ] - C(G)MM^) and a 
possible dependence on the unphysical mass of the e scalar. As emphasized in [35], one 
should properly allow the e scalars to have masses and mass counterterms in the component 
field approach to DRED, and the result given there is correct. However, one can treat the 
subtractions in DR in such a way that the two-loop (3 functions do not depend on the un- 
physical e scalar mass, and so that the expressions relating the pole masses to the running 
masses of the scalars also do not depend on the e scalar masses. (This is related to the 
prescription used in [35] by a simple coupling constant redefinition for the scalar masses, 
of the type mentioned there.) The corrected results we have given here correspond to this 
prescription, and we are now in agreement with the authors of [34] and [35], to whom we 
are grateful for consultations. The subtleties involved will be reported on elsewhere. The 
results in the special case of the MSSM have also been corrected accordingly. 
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